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Abstract 

We construct a new model of the quantum oscillator, whose energy 
spectrum is equally-spaced and lower-bound, whereas the spectra of 
position and of momentum are a denumerable non-degenerate set of 
points in [—1, 1] that depends on the deformation parameter q G (0, 1). 
We provide its explicit wavefunctions, both in position and momentum 
representations, in terms of the discrete q-Hermite polynomials. We 
build a Hilbert space with a unique measure, where an analogue of 
the fractional Fourier transform is defined in order to govern the time 
evolution of this discrete oscillator. In the limit when q ^ \~ one 
recovers the ordinary quantum harmonic oscillator. 
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1 Introduction 



Several algebraic constructions have been proposed in the hterature to de- 
scribe various extensions of the quantum harmonic oscillator. These con- 
structions are based on various deformations of the standard oscillator Lie 
algebra, or different associative algebras. In most of these models it is dif- 
ficult to construct a theory for such oscillators, which is as complete as the 
well-known treatment of the standard harmonic oscillator in quantum me- 
chanics. Namely: a canonical complementarity between position and mo- 
mentum, explicit forms for the wavef unctions, and a coherent description of 
time evolution. 

The earliest model, generally called the g-oscillator, was proposed by 
Macfarlane [1] and Biedenharn [2] on the basis of raising and lowering eigen- 
states of a Hamiltonian with the g-deformed commutator Og — qaqa'^ . A 
theory of this oscillator has been elaborated (see, e.g. [3]-[6]); yet, it has not 
been clear how to construct position and momentum operators satisfying the 
basic commutation relations with a Hamiltonian to characterize infinitesimal 
harmonic motion. This may be one of the reasons why this q'-oscillator has 
not proven attractive for many physicists. 

The postulates we use to define oscillator models are the following [7, 8]: 

1. There exists an essentially self-adjoint position operator Q, whose spec- 
trum X is the set of positions {x} of the system. 

2. There exists a self-adjoint and compact Hamiltonian operator H, whose 
commutator with position defines the momentum operator P, 

[i/,g]=:-iP, (1) 

and corresponds to the first Hamilton equation (i = a/— T); the commutator 
of the Hamiltonian with momentum returns the position operator 

[H,P] = iQ, (2) 

that corresponds with the second Hamilton equation, and which characterizes 
the oscillator dynamics. Equivalent to (l)-(2), one can propose the Newton- 
Lie equation as [iJ, [H, Q\\ = Q. The set of momentum values of the system 
is the spectrum of P, which is equal to that of Q because (l)-(2) generate a 
rotation between these two operators (to be written below). 
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3. The three operators, Q, P and i7, close into an associative algebra^ i.e., 
they satisfy the Jacobi identity, 

[P, [H, Q]] + [Q, [P, H]] + [H, [Q, P]] = 0. (3) 

We note that the basic commutator [Q, P] has not been defined. Due to 

(I) -(2), the only restriction imposed by the associativity condition (3), since 
the first two summands are identically zero, is that [Q, P] must commute 
with H, and thus be constant under the oscillator motion. This indicates 
that each distinct choice of the basic commutator [Q, P] will yield a dis- 
tinct model for the oscillator. If the choice is the Heisenberg commutator 
[Q,P] = hi, one has the standard four-generator oscillator Lie algebra Ha ~ 
span {H, Q, P, 1} of quantum mechanics (containing the Heisenberg algebra 
Tis = span {Q,P,1}). In one previous work [9], with the purpose of endowing 
Q with a finite set of position eigenvalues x G {—j, ■ ■ ■ ,j} (we write 
xI'Lj), the basic commutator was taken to be [Q, P] = i[H — {j + |)1] =: i J3, 
in a matrix representation of the Lie algebra so(3) = su(2) = span {Q, P, J3} 
of spin j, so the three operators have the same finite spectrum x\^_j. This 
model, called the su(2) oscillator, has been applied to study the parallel pro- 
cessing of finite signals and pixellated images [10]. The general conditions to 
include oscillator dynamics in associative algebras were given in [7]. 

In a previous paper of the same authors [8] , a g-algebraic associative struc- 
ture was proposed on the basis of the quantum algebra sug(2), the Hamilto- 
nian having a lower-bound equally-spaced spectrum. Using a non-standard 
basis to define position and momentum operators that allowed analytic ex- 
pressions, their spectra was determined to be a finite set of non-equally spaced 
points Xs = sinh(sfi;)/ sinh ^k, with s\-'_j and q = e"" G (0, 1). Also, explicit 
expressions were obtained for the wavefunctions in position and momentum, 
in terms of the dual g'-Kravchuk polynomials, related by a fractional finite 
Fourier-g-Kravchuk matrix transform, and a natural representation in a sui 
generis phase space. The present paper is a continuation of the research in 

[II] that constructed quantum oscillators with continuous bounded spectra 
for the position and momentum operators. 

In this paper we build an oscillator model on the basis of the Fock repre- 
sentation of a quadratic associative algebra which is a g-deformation of the 
standard oscillator algebra; we denote this by 'DTii, which will be defined 
in Section 2, while the physical interpretation of the participant operators 
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that characterizes this model are set forth in Section 3. The spectrum of 
the Hamiltonian in the algebra is lower-bound and equally spaced, as in its 
standard counterpart. The position spectrum and wavef unctions, orthonor- 
mal under a specific scalar product over positions, are obtained explicitly 
in Section 4, and the momentum wavefunctions in Section 5. This model 
can be characterized for having a space of positions given by an infinite non- 
degenerate point set contained in the interval [—1, 1]. The coordinate and 
momentum realizations of the oscillator are given in Sections 6 and 7; in 
terms of these we can write the harmonic oscillator motion in Section 8, with 
a summation kernel is a new, which is a (denumerable infinite) relative of 
the fractional Fourier transform of the standard case. Concluding remarks 
are offered in Section 9. 

We use the notations that are standard in the theory of basic hypergeo- 
metric functions and g-orthogonal polynomials (see, for example, [12]), and 
we assume throughout that g is a fixed real number in (0, 1). 

2 The quadratic algebra Vq Ha 

We define the algebra VgH^ as the associative algebra generated by a vector 
basis of elements /+, Jq, satisfying the following commutation relations: 

[lo, /±] = ±/±, [/+, /-] = - (1 + g) q''° . (4) 

Equivalently, introducing Ii :— 7+ -|- 7_ and I2 '■— i(7+ — 7-), we can charac- 
terize this algebra by 

[7o, h] = -172, [lo, h] = i/i, [h, I2] = -2i{q'° - {1 + q) q''') • (5) 
The first relation in (4) can be written in the form 

q'°I±q-'°^q^'l±. (6) 

This relation and the fact that both q^° and q^^" appear in the second rela- 
tion of (4) shows that VgH^ — span {7i, 72, g^^"} is a quadratic associa- 
tive algebra. This is a g-deformation of the oscillator algebra because 
lim^^i- VqHi — Hi. Indeed, in the limit lim^_»i- we obtain from (5) the 
relations 

[7o, h] = -172, [lo, I2] = ill, [h, h] = 2i, (7) 
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which are equivalent to the defining relations of 7^4. 

We are interested in the Fock representation of the algebra VgH^; this 
is an irreducible representation constructed on a Hilbert space with the or- 
thonormal basis of vectors e„, n G {0, 1, 2, ■ • • } {i.e., In this represen- 

tation, using the 'box' g-number [a]q := (1— g), the operators of the 
algebra act by raising and lowering the number n of e„, 

/+e„ = Y^g"+i[n+l]ge„+i, /_e„ = y^g"[n]g e„_i, (8) 

/oe„ = ne„, i.e., q^"en = q''en, (9) 

and the hermiticity conditions = /_ and Jq = /q are satisfied. 

In order to have a functional realization of this representation, we consider 
the space V of all polynomials in one supplementary variable y, and introduce 
its basis of monomials 

^n(n-l)/4 

Cn ^ e„(y) := c„ y", c„ = —j^, n\^, (10) 

q)n 

where (a; q)n '■— (1 — a)(l — aq) ... (1 — ag"~^) and (a; q)o — 1. Acting on 
analytic functions f{y) G V, the Fock representation can be written in terms 
of the scale Ta and g-difference Dg operators, 



q'°^Tg, I_= [q{l-q)]'/'Dg; (11) 

Tam=f{ay), Dgf{y) = M^^^. (12) 

This realization of the algebra is equivalent to that in (8) -(9), with the 
functions e„(y) playing the role of the basis elements e„ as eigenfunctions of 

the weight operator Jq. 

Let us now introduce a scalar product into the space of polynomials V. 
This scalar product is of a Fisher-type scalar product and is given by the 
formula 

{h{y)j2{y))^fi{Dg)fl{y)\y=o, 

where Dg := (1 — q)Tg-iDg, fi and /2 are polynomials, and fl denotes the 
polynomial /2, whose coefficients are replaced by their complex conjugate 
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ones. In this formula, we have the action of the difference operator upon the 
polynomial f^- Then it is directly verified that 

{en: en') = Sn,n': ri^u'l^. (13) 

Closing the space V with respect to this scalar product we obtain a Hilbert 
space that we denote by H. The space IH consists of functions 

oo oo oo 

n=0 n=0 n=0 

where a„ = &nC„, and c„ are determined by (10). Since (e„, e„/) = 5n,n' by 
definition, for /i(y) = J2n=o (^ny"" and /2(y) = X)^o have 

oo , 
n=0 I "I 

This means that the Hilbert space H consists of analytic functions /i(y) = 
J2n=o ^nV"^ such that 

< oo. (16) 



3 Assignment of observables to generators 

The discrete oscillator is a class of oscillator models that depend on the 
parameter q G (0, 1), and based on the Fock irreducible representation of the 
algebra T>qH4, where the physical observables are assigned to the spectra of 
self- adjoint generators of the algebra in the following way: 

position: Q := ^(l^^/i, (17) 
momentum: P :— \/ {l—q)/ql2, (18) 
Hamiltonian: H := Iq + |i. (19) 

Then, due to (5), H exhibits the Hamiltonian oscillator commutation rela- 
tions (l)-(2) with Q and P, and determines the basic commutator [Q,P], 
through 

[H,Q] = -iP, [H,P] = iQ, (20) 
[Q,P] = 2i{l-q-')[q"-'/'-{l + q)q'"-'] =:iF{H). (21) 



6 



The operator F{H) defined in (21) commutes with the Hamiltonian H and 
is therefore also diagonal in the Fock basis {e„},^Q in (10), 

F{H) e„ = 2(1 - q-^) [q^ - (1 + g) q^^] e„. (22) 

This basis of H thus consists of eigenfunctions of a Hamiltonian with equally- 
spaced eigenvalues, 

i7e„ = (n+|)e„, n|~, (23) 

coinciding with the energy spectrum of the standard quantum harmonic os- 
cillator. 

Prom (20), the time evolution of the discrete oscillator position and mo- 
mentum operators is produced by 

exp(iri/) = e'^/^ exp(iT7o), (24) 

and results in the harmonic motion 

(«J:|)^e-(«)e-^(_- -)(^). (25) 

which for r G [0, 2ti) forms a group U(l) of inner automorphisms of the 
pair of operators Q and P, that we interpret as rotations of a phase plane 
around its origin (still to be studied for this model, see [8]). The phase e'"^/^ 
is due to the energy ^ of the ground state cq, while exp(ir/o) is the discrete 
counterpart of the fractional Fourier transform for this model, to be seen in 
Section 8 below. 

The associative algebra VgH^ is a g-deformation of the standard Heisen- 
berg-Weyl algebra, where I± in (8) are recognized as the raising and lowering 
operators 

a+:=/+ = iyi^(g-iP), a,:=/_ = iyi^(g + iP). (26) 

From (4) it then follows that 

[ lim Oq, lim a+] = lim Ul+q) q^^° - q^"] = 1, (27) 

so we recover the standard oscillator. This is the place to emphasize the 
difference between the Macfarlane-Biedenharn g-oscillator [l]-[6], which is 
defined in terms of g-commutators aqO^ — q ttq, and our discrete oscillator, 
which is formulated exclusively with ordinary commutators. 
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4 Spectrum and eigenfunctions of the posi- 
tion operator 



A direct calculation using (8) shows that in the Fock eigenbasis {e„}^Q of 
the Hamiltonian H, the position operator Q — q"^ — lli acts as 

Qen = - e„+i + en-i • (28) 

Since — q^^'^)\ < 1 for the norm ||Q|| of Q docs not exceed 1, 

and hence it is a bounded operator whose eigenvalues {x} will lie in the real 
interval [—1,1]. 

4.1 Eigenvectors of position 

To find the eigenvectors ipxiu) and the spectrum {x} = A" of the position 
operator Q, 

QMv) = xi/j^iy), xeX, (29) 
we represent ipx{y) in the form of a linear combination of the monomials (10), 

oo 

tpxiy) = ^Pn{x)en{y), (30) 

n=0 

where Pn{x) are coefficients depending on the points of the spectrum x & X. 
When we substitute the expansion (30) into the equation (29), we obtain 

oo 

J2 Pn{x) Vg"(l - ?"+!) e„+i + pn{x) V?"-i(l-g") e„_i 

"=° oo (31) 

^x^Pri{x) e„. 

n=0 

From here we find the following three-term recurrence relation for the coef- 
ficients Pn{x) in (30), 

xpn{x) = y/q^{l-q^+^)pn+i{x) + y/q^-^{l-q^) Pn-i{x) , (32) 
starting with p-i{x) = and po{x) := 1 setting the common constant factor. 
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We see from (32) that the coefficients Pn{x) in (30) are polynomials in x of 
degree n, which can be evaluated uniquely. To solve the recurrence relation, 
we make the substitution 

Pn{x) = {q;q)-'/\-^--'^/%{x), (33) 

which turns (32) into 

xp„(x) =p„+i(x) + ?"-'( 1 - q")pn-iix). (34) 

Comparing this with the recurrence relation for the discrete g-Hermite poly- 
nomials of type I, given by 20i basic hypergeometric polynomials [13, Eq. 
(3.28.3)], 

hn{z;q) := g'^("-^)/Vi(g-", ^-^0; -gz), (35) 
zhn{z;q) = hn+i{z;q) + q''-\l-q'')hn-i{z;q), (36) 

we establish that Pn{x) — hn{x; q). We can thus write the coefficient polyno- 
mials in (30) as 

Pnix) = iq;q)-'/'q-^-'yXXx;q)- (37) 
From (32) follows that these polynomials have definite parity: Pn{x) = 

i-lTPnix). 

Collecting these results we write the eigenfunctions ip^ iu) of the position 
operator Q as 

oo 

UV) = J2i^;q)n'^'q-''^''-'^^'hn{x;q)eM 

n=0 

oo 

= ^{Q;Q)n^hn{x;q)y'', 

n=0 

{xy;q)oc ' 

where in the last expression we use the symbol (a; q)oo '■= Y[^=o (-'- ~ '^9'") 
the summation formula in [13, Eq. (3.28.11)]. Because of the convergence of 
{y^'i<f)oo in (40) for the basis {'il^x{y)}xex, we must restrict the domain of 
definition of functions f{y) e H to the open disk \y\ < 1. Then the condition 
(xy; q)oo < 1 is fulfilled automatically since we saw that the eigenvalues x & X 
of Q are contained in the interval [—1, 1]. 
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(38) 
(39) 
(40) 



4.2 The spectrum of position 

The spectrum of the self-adjoint position operator Q ~ Ji = /+ + /_ can be 
found from the series (38); from (28) we see that in the basis {e„(|/)} 
operator Q is a self-adjoint Jacobi tridiagonal matrix of the form 



oo j-Up, 



Q 



/bo 

ao 



V : 



ao 

h 







02 












03 



\ 



o„ 7^ 0. 



(41) 



We can now use the theory of these matrices from [14, Chap. VII], (see also 
[6]) to connect their spectra with the corresponding measures for orthogo- 
nal polynomials. In this vein, we note that in the Fock basis the position 
eigenfunctions ipxiv) are expanded in terms of the basis elements {e„(i/)}^ 
with the polynomial coefficients Pn{x) in (30), which are given in terms of 
discrete gr-Hermite polynomials of type I in (37). According to the results 
in [14, Chap. VII] , these polynomials are then orthogonal with respect to a 
spectral measure dfi{x) of the operator, which is unique up to a constant 
factor, on a set A" C 3? that is the simple spectrum of Q. 

In finding the spectrum of the position operator Q, we recall that the 
discrete g-Hermite polynomials hn{x; q) obey the orthogonality relation 



q X ]q)oo hk(x; q) hm(x; q) dqX 



-- 4,m (i-'?)(9^; ?^)c 

CO 



-l;g)oo {q; q)mq" 



m(m-l)/2 

1 

-q;q)liq;q)mq"'^"'-'^/', 



(42) 



where J^^ f(x)dqX is the symbol of the g-integral (see [13, Eq. (3.28.2)]). 
This orthogonality relation can be written in the form of a sum [12], 



n=0 



m(m-l)/2 



= 25fe,„ (g; q)^ (-g; q)^ (g; q)^ g" 
This means that the spectrum A" of Q is the simple set of points 

A' = {g",-g"; n\^}, 



(43) 



(44) 
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and that the corresponding eigenfunctions are 

i/jqniy), i/j-qniv), (45) 

given by (38)-(40). The spectrum of Q is discrete, which means that the 
eigenfunctions ip±qn{y) form a denumerable orthogonal basis in the Hilbert 
space H; we note that X C [—1, 1] has a unique accumulation point that 
does not belong to the set. 



4.3 Normalization of the eigenfunctions 

The eigenfunctions of Q were determined only up to constant factors, so we 
proceed to normalize the eigenfunctions {i^±q''{y)}'^ in their form (39). Prom 
(38) and the orthogonality of the basis {e„(y)}^Q we obtain 

oo^ ^-n(n-l)/2 

{i^x{y),ipx'{y))m = 2^ —7 ^ — Kix; q) hn{x'] q), (46) 

n=0 

where x and x' take values in X = {q^, —q^; s\^}. We can calculate this 
sum as follows: we build the functions 



We can see the hn{^q^', q) as the elements (n, ±s) of a matrix of numbers 
(integer n > numbering rows and s > numbering columns), written as 



R(g^;?)}~=o 



(49) 



Columns of this matrix are orthonormal due to the orthogonality relation 
(43) for the discrete q'-Hermite polynomials hk{z;q). In the infinite dimen- 
sional case, the orthonormality of columns does not immediately lead to the 
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orthonormality of rows. But in accordance with the reasoning of Ref. [15], 
one can state that rows of this matrix are also orthogonal, i.e., 



n=0 

oo 

n=0 

oo 

J2'^n{-q']q)hn{-q"' 

n=0 

Substituting (47)-(48) into (50)-(52), we obtain 

oo 

E 



q) = Ss,s', 

q) = 0, 

q) = Ss,s'- 



(50) 
(51) 
(52) 



{q''+';q'Uq' ^ K{±q^ ; q) K{±q^' ; q) 



2(g; g)oo (-g; q)l ^ g"^"-^)/^ 



(53) 



where one has to take only the upper or only the lower signs. Returning to 
the scalar product in (46), we find 



/ / / N / / N\ r 2(g;g)oo (-g;g)L 



(54) 



We thus arrive at the functions 



•,2s+2. „2 



2(?; g oo(-?; ? So 



which are orthonormal under the scalar product (46) in H, 
(*^(y),*a;/(y))H = x,x'eX. 



(55) 



(56) 



5 Spectrum and eigenfunctions of the mo- 
mentum operator 

The momentum operator P I2 — i(/+ — /_) acts on the basis {e„(y)}2° as 
Pe„ = i(Vg-(l-r+i)e,+i - ^g-i(l-gn) (57) 
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[cf. (28)]. When we change this basis to another {e„}[^ with e„ = i^e^, one 
can see that in the new basis the momentum operator P acts as a matrix 
with the same coefficient elements as the position operator in Section 4 on 
the former basis. This means that the spectrum of momentum P coincides 
with the spectrum of position Q, namely, SpecP = X, where X is given in 
(44). Similarly, the cigcnfunctions of momentum P can be found in the same 
way as the eigenf unctions of Q, by using the basis {e„}Jf. 

Let (f)p{y) satisfy P(f)p{y) — p(f)p{y), an eigenfunction of P corresponding 
to the eigenvalue p, with an expansion in the mode eigenbasis {e„(y)}^ given 

by 

oo 

My) ^^9n(p) en(y), (58) 

n=0 

where Qnip) are coefficients depending on the momentum p ^ X. Repeat- 
ing the process of the previous section, one derives a three-term recurrence 
relation for the polynomials Qnip) and concludes that 

( \ -n ( \ i''hn{p;q) . s 

9n[p) = 1 Pnip) = ; ,1/2 , (59) 

where hn{z\ q) are the discrete g-Hermite polynomials of type I from Section 
4. Hence, the eigenfunctions of momentum are 

= '7/f?\y, 'r.{y) (60) 

pe A' = {g^-g- s\^}. (62) 



n=0 

{y^;q)c 
(iyp; q) 



To find the last two expressions we have used the same method as in the case 
of eigenfunctions of position in (38)-(40). 
The normalized eigenfunctions of P are 



'oo 



satisfying {^^{y),^^,{y))M = 4,a;', x,x' e X. 
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6 Coordinate realization of the discrete oscil- 
lator 

In Section 3 we constructed a realization of the discrete oscillator on the space 
of analytic functions in the supplementary variable y with the assignment 
(10). It is natural to look for a realization of the oscillator on the space of 
functions in the position coordinate x E X . 

Let L'^{X) be the Hilbert space of square-summable functions over x E X 
(the set of positions of the discrete oscillator) , with the scalar product 

^ oo 

(r;r)oo(-i;g)oo (64) 



X (/i(g")/2*(?'^) + /i(-?'^)/2*(-?")), 



where * stands for complex conjugation. 

Since the discrete g-Hermite polynomials are associated with the deter- 
minate moment problem (see, for example, [6] for the description of this 
association), the set of polynomials {Pn{x)}'^ in (37) constitute a complete 
set of orthonormal functions in the Hilbert space L'^{X). 

We construct a one-to-one hnear isometry Q from the Hilbert space H, 
onto the Hilbert space L'^{X), given by 

Q : e 9 eiy) ^ f{x) = {e{y),My))M e L\X), (65) 

where ipxiu) are eigenfunctions (40) of Q. It follows from (38) that 

H 3 e„(y) (e„(y), '0x(z/))h = Pn{x). (66) 

That is, Q maps the basis {e„(y)} of H, which is orthonormal under the 
scalar product (54), onto the basis {pn{x)} of L^(A'), which is orthonormal 
under (64); this means that $7 is a one-to-one isometry. 
In L^(A'), the operator Q acts through multiplication, 

Q fix) ^x fix). (67) 

Indeed, since Qipxiu) = xip^iy) for Jle(?/) = fix) = (e(y), V'j;(z/))h, we have 

n-.Qeiy) Q fix) = {Qe{y),i)x{y))m ^^^^ 
= {eiy),Qi)xiy))m = {eiy),xijjxiy))u = xfix). 
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We can find the action of Q, P, and H on the basis elements {p„(a;)}^Q 
of the Hilbert space L^(A'). According to the recurrence relation (32), which 
follows from the recurrence relation for the discrete q'-Hermite polynomials 
hn{z; q), we have for the position operator Q that 

QPn{x) = + y/q^-\l - q'^) Pn-l{x) . (69) 

It follows from formulas (3.28.7) and (3.28.8) in [13] that the momentum 
operator P acts on the Hilbert space L'^{'^) through 

P = -i(l - ,) (b, + -.j^^D,-, (,V; , (70) 

where (g^x^; g^)oo is the multiplier in the orthogonality measure in the scalar 
product (64). In particular, P acts on the basis functions {Pn{x)}'^^Q as 

Pp^{x) - i^yq^{l-q^+^)pn+i{x) - iVg"-Hl-?")Pn-i(^)- (71) 

Finally, the Hamiltonian H acts on the basis polynomials Pn{x) of the 
Hilbert space L'^{X) as 

Hpn(x)^(n+l)pn(x). (72) 

Indeed, according to (23) and (66) we have 

Hpn{x) = {Hen{y),^x{y))m ^^^^ 
= {n+ l){en{y),i'x{y))w = (n + i)p„(x). 

It is interesting to see the lowest discrete oscillator modes Pn{x) in (37) 
and (72), both as continuous functions of x G (0, 1] and their values at the 
orthogonality set X, for various values of q. While the continuous functions 
exhibit strong oscillations (increasing with n and 1/q), their values on X{q) 
remain well bounded. As q ^ 1~ , the corresponding points in X densify, 
evincing the resemblance of the discrete oscillator with the standard oscillator 
wavefunctions. This form of convergence should be studied further, since a 
change of scale appears necessary as well as a discrete measure that becomes a 
continuous Riemann integral in the limit. One analogue for this limit appears 
in [16, Fig. 4], where the Meixner functions that describe the discrete model 
converge to the Laguerre- Gauss modes of a radial oscillator; in that case 
though, the limit is from hyperboloids to the cone in the three-dimensional 
space of the Lie algebra su(l, 1), and is not a g-deformation. 
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7 Momentum realization of the discrete os- 
cillator 



Consider the Hilbert space Li^{V) of square-integrable functions f{p) in the 
momentum coordinate p in the oscillator with the same scalar product as in 
(64), where V = X is the spectrum of the momentum operator P, coinciding 
with the spectrum of Q. The coefficient polynomials gn{x) in formula (59) for 
the eigenf unctions of momentum, (j)p{y) in (58), constitute an orthonormal 
basis in L^iV). 

To formalize this consideration, we construct, as in the previous Section, 
a one-to-one linear isometry Q from the Hilbert space H onto the Hilbert 
space L'^{V), given by 

Q: m3 e{y) ^ f{p) := (e(y), My)h e L'{V), (74) 

where 4>p{y) are the eigenfunctions of momentum P in (60)-(62). [Compare 
with (65) requiring the position eigenfunctions ipx{y) in (38)-(40).] Prom 
here it is evident that 

H 3 Cniy) (eniy), 0p(|/))h = gn{p): (75) 

that is, is a one-to-one isometry and maps the orthonormal basis {e„(y)}^ G 
H onto the orthonormal basis {gn{p)}'o' G L^i'P)- 

The momentum operator P acts on LF'iV) as a multiplication operator 
on all functions of p, 

Pg{p)^pg{p). (76) 

The action of Q, -P, and H on the basis of polynomials gn{p) can be found 
in the form of the recurrence relations 

Qgn{p) = V g"(l - g"+^) g n+iip) + v ^-^l - g") g n-i(p), (77) 
Pgn{p) = iVg"(l - g"+i) ^„+i(p) - iV?--i(l - (t) 9n~i(p), (78) 
Hgnip) = (ri+|)^„(p). (79) 

8 Harmonic evolution in position space 

According to (24) and (25), the action of the operator exp(iriJ) is the time 
evolution of the discrete oscillator. On the basis (10) of functions {e{y)}^ 
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that is orthonormal with respect to the scalar product (13), this action is 

e-^ e„(t/) = e-/2 e^^ e„(t/) = e^^^+^Z^)- e„(y). (80) 

The operator exp(iri7) also acts on the Hilbert space L'^{X), which is char- 
acterized by the scalar product (64). Now consider the isometry between 
these two spaces, 

H 9 e{y) ^ fix) (e(t/), My))^ G L2(A'), (81) 

that maps functions e{y) onto functions f{x) of the discrete position coordi- 
nate X G {—q^,q^}'^ = X. Then to exp{iTH)e{y) G H there corresponds a 
function exp(irif)/(±g*) of a; = ±g*, 

e-^/(±g^) = (e-^e(y),V±,^(y))H=(e(2/),e-^V'±,.(|/))H (82) 

oo 

= 5Z(e(l/),en)H(en,e~'^-^V'±g.(l/))H (83) 

n=0 

oo 

= X]^^^^)'^")h^^'"'^^"'^±9^^1/))h (84) 

n=0 

oo CO 

n=0 m=0 

+ {e{y), ^-,^{y))M{^-qm{y), e„)H) (85) 

oo 

= J] (ir^(± g^ g™)/(D + K^{± q% -g")/(-g™)) ;(86) 

m=0 

where, according to (38), 

oo 

X-(±g^±g-) = c,X(V'±,^(y),e„)He-("+V2)(e„,V'±,.(y))H (87) 



n=0 



°° ^n(n-l)/2 

c.e^^/^ E / A ^"(± e'"' ^n(± ?) (88) 



and 



oo 



2(g;g)oo(-?;?) 
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2 ' 

oo 



(89) 



Being elements of L'^{X), the functions f{x), x G Af, enter into the scalar 
product (64) as a sum with the weight function (a;^g^; g^)oo- Because we 
intend to interpret these as wavef unctions of a quantum oscillator, and dis- 
play them for comparison with their standard shapes, we should absorb this 
weight into the functions. 

Fix) = VW-,Q')oof(x), (90) 

so that their scalar product acquires, from (64), the 'standard' form, 

(r;r)oo(-i;g)oo ^ 

X = x{±,n) = ±q'^ e X, n\^. (91) 

The matrix elements of operators will be correspondingly rescaled by (90). 

The oscillator evolution (82)-(86) can be used to define the fractional 
discrete Fourier transform on the functions F{x), x E X, rescaled as in (91). 
We note that the fractional Fourier integral transform for angle r differs from 
the standard harmonic oscillator evolution by a phase e'"^/^ that is due to the 
ground energy | in the oscillator, so the fractional Fourier transform is 

$(r) := e-^^/^exp(iT//). (92) 

Its action on the rescaled discrete wavefunctions will have the form 

$(r)F(a;) = ^^{x,x' ^t) F{x'), (93) 
x'ex 

where x' G and K'^{x,x') is given by (87)-(88) for x = ±q'^ and 
x' = ±g™. Unfortunately, the bilinear generating function (88) of discrete 
qf-Hermite polynomials of type I could not be summed to a closed form. 

9 Concluding remarks 

We constructed a model of the harmonic oscillator that can be realized on 
bases of coordinate and momentum Hilbert spaces, and its energy modes ex- 
pressed in terms of discrete g-Hermite polynomials of type I. The spectrum 
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of the Hamiltonian coincides with that of the standard harmonic oscillator 
in quantum mechanics, while the position and momentum operators in this 
model have discrete, denumerably infinite spectra that depend on the exten- 
sion parameter q contained in the interval [—1, 1]. 

Contrary to other models that use discrete g-Hermite polynomials [17, 18, 
19], our models (the present one and that in Rcf. [8]) fulfill the basic Hamil- 
ton equations in the form [H, Q] = —iP and [H, P] = iQ, with standard 
commutators — and not g-commutators [1, 2]. Because of this important 
circumstance, the time evolution of the model is a Lie group, which but for 
a phase is that of fractional Fourier transforms associated with this model 
[20] . This discrete oscillator is a new and non-trivial deformation of the stan- 
dard quantum harmonic oscillator; it allows the extension of other standard 
concepts of phase space, such as coherent states [8], that will be examined 
elsewhere. 

We believe that the discrete oscillator model can appropriately describe 
discrete quantum systems on bounded point lattices, and also contribute 
significantly to the general theory of special functions. 
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